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We calculate the form factors for weak decays of Br s -j and Dr s \ mesons to light pseudoscalar and 
vector mesons. To reveal the intimate connection between different decay modes and to be able to 
perform the calculations in the full physical q 2 -region we use a relativistic dispersion approach based 
on the constituent quark picture. This approach gives the form factors as relativistic double spectral 
representations in terms of the wave functions of the initial and final mesons. The form factors have 
the correct analytic properties and satisfy general requirements of nonperturbative QCD in the 
heavy quark limit. 

The disadvantages of quark models related to ill-defined effective quark masses and not precisely 
known meson wave functions are reduced by fitting the quark model parameters to lattice QCD 
results for the B — > p transition form factors at large momentum transfers and to the measured total 
D — ► (K,K*)lv decay rates. This allows us to predict numerous form factors for all kinematically 
accessible q 2 values. 

PACS numbers: 12.39.Ki, 13.20.He, 13.20.Fc 



I. INTRODUCTION 



The knowledge of the weak transition form factors of heavy mesons is crucial for a proper extraction of the quark 
mixing parameters, for the analysis of non-leptonic decays and CP violating effects and, related to it, for a search of 
New Physics. 

Theoretical approaches for calculating these form factors are quark models QCD sum rules [p|— p~2| , and lattice 

QCD g|-|l|. Although in recent years considerable progress has been made, the theoretical uncertainties are still 
uncomfortably large. An accuracy better than 15% has not been attained. Moreover each of the above methods has 
only a limited range of applicability, namely: 

QCD sum rules are suitable for describing the low q 2 region of the form factors. The higher q 2 region is hard to get 
and higher order calculations are not likely to give real progress because of the appearance of many new parameters. 
The accuracy of the method cannot be arbitrarily improved because of the necessity to isolate the contribution of the 
states of interest from others. 

Lattice QCD gives good results for the high q 2 region. But because of the many numerical extrapolations involved 
this method does not provide for a full picture of the form factors and for the relations between the various decay 
channels. 

Quark models do provide such relations and give the form factors in the full q 2 range. However, quark models are 
not closely related to the QCD Lagrangian (or at least this relationship is not well understood yet) and therefore have 
input parameters which are not directly measurable and may not be of fundamental significance. 

In this situation it becomes evident that a combination of various methods will be fruitful. It is the purpose of this 
paper to obtain this way reliable predictions for many decay form factors in their full q 2 ranges. 

To achieve this goal, one needs a general frame for the description of the large variety of processes. This can be only 
a suitable quark model, because only a quark model can connect different processes through the soft wave functions of 
the mesons and describes the full q 2 range of the form factors Q The predictions of this model will be much improved 
by incorporating the results from the more fundamental QCD based methods. 

Historically, constituent quark models have been first to analyse the meson transition form factors. Although 
a rigorous derivation of this approach as an effective theory of QCD in the nonperturbative regime has not been 
obtained, relativistic quark models work surprisingly well for the description of the meson spectra and form factors. 
Moreover, constituent quark models provide so far the only operative method for dealing with excited states. 



1 One of the first steps in combining various approaches in order to obtain predictions for the form factors for all q 2 has been 
done in |l^| where a simple lattice-constrained parametrization based on the constituent quark picture of Ref. jl^] and pole 
dominance has been proposed. 
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1. The physical picture 



Constituent quark models are based on the following phenomena expected from QCD: 

i) chiral symmetry breaking in the low-energy region provides for the masses of the constituent quarks; 

ii) the strong peaking of the soft (nonperturbative) hadronic wave functions in terms of the quark momenta with 
a width of the order of the confinement scale; and 

hi) the dominance of Fock states with the minimal number of constituents, i.e. a qq composition of mesons. 

An important shortcoming of previous quark model predictions was a strong dependence of the results on the special 
form of the quark model and on the parameter values. 

The goal of this paper is to demonstrate that once (a) a proper relativistic formalism is used for the description 
of the transition form factors and (b) the numerical parameters of the model are chosen properly (we discuss criteria 
for such a proper choice below), the quark model yields results in full agreement with existing experimental data and 
in accord with the predictions of more fundamental theoretical approaches. In addition, our quark model allows to 
predict many other form factors which have not yet been measured. 

2. The formalism 

For the description of the transition form factors in their full q 2 range and for various initial and final mesons, a 
fully relativistic treatment is necessary. We therefore choose a formulation of the quark model which is based on the 
relativistic dispersion approach [ ft8| and thus guarantees the correct spectral and analytic properties. 

Within this model, the transition form factors are given by relativistic double spectral representations through 
the wave functions of the initial and final mesons both in the scattering and the decay regions. These spectral 
representations obey rigorous constraints from QCD on the structure of the long-distance corrections in the heavy 
quark limit. Namely, the form factors of the dispersion quark model have the correct heavy-quark expansion at 
leading and next-to-leading 1 / toq orders in accordance with QCD for transitions between heavy quarks [fL9|j2(]] . For 
the heavy-to- light transition the dispersion quark model satisfies the relations between the form factors of vector, 
axial- vector, and tensor currents valid at small recoil |2l[. In the limit of the heavy-to-light transitions at small q 2 
the form factors obey the lowest order 1 /ttiq and 1 /E relations of the Large Energy Effective Theory |||] and provide 
the pattern of higher-order symmetry-violating effects. 

Another important advantage of the dispersion formulation of the quark model is the fact that one directly obtains 
the form factors in the physical decay region q 2 > 0. No numerical extrapolation from space-like values is required. 

3. Parameters of the model 

In previous applications of quark models the transition form factors turned out to be sensitive to the numerical 
parameters, such as the quark masses and the slopes of the meson wave functions. As proposed in Ref . ]25} ] , the way 
to control these parameters is to use the results of lattice calculations at large q 2 as 'experimental' inputs. In (24] the 
b and u constituent quark masses and slope parameters of the B, 7r, and p wave functions have been obtained through 
this procedure. 

We now consider in addition the charm and strange mesons. To determine the slope parameters for the charm and 
strange meson wave functions and the effective mass values m c and m s we use here as input the measured total rates 
for the decays D — > (K, K*)lv. By fixing the parameters in this way we overcome important uncertainties inherent in 
quark model calculations. Indeed, with these few inputs we can give numerous predictions for the form factors for the 
transitions of the heavy and strange heavy mesons D, D s , B, and B s into light mesons which nicely agree at places 
where data are available. 

The paper is organized as follows. Section II briefly presents the main features of the double spectral representations 
of the form factors. In Section III we determine the numerical parameters of the model and give the predictions of 
the form factors. Section 4 contains our conclusions. 



II. MESON TRANSITION FORM FACTORS 

The long-distance contribution to meson decays is contained in the relativistic invariant transition form factors 
of the vector, axial-vector and tensor currents. The amplitudes of the M\ — ► M.% transition induced by the weak 
Ii —> Qi quark transition through the vector = 917^92, axial-vector = 9i7M7 5 g2, tensor T^ v — qi<J IJ , I/ q2, and 
pseudotensor T^ v = q\cr iiv ^q2 currents, have the following covariant structure [ pl| 
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< P(M 2 ,p 2 )\V fl (0)\P(M uPl ) > = f+(q 2 )P^ + f-(q 2 )q^ 



< V(M 2 ,p 2 ,e)\V^0)\P(M llPl ) > = 2g(q 2 )e^ a0 e w p? pf, 

< V(M 2 ,p 2) e)\A fl (0)\P(M 1 ,p 1 ) > = ie* a [f{q 2 )g^+a + {q 2 ) Pla P^a^{q 2 ) Pla q^ 
< P{M 2 ,p 2 )\T llv {Q)\P{M 1 ,p 1 ) > = ~2is{q 2 ) (p lM p 2 „ nu p 2fl ), 

< V(M 2lP2 ,e)\T^(0)\P(M llPl ) > = ie* a [ 9+ {q 2 )e^ a pP^ + g-(q 2 )e^ a ^ + g a {q 2 ) Pla ^ v ^p{p\ }, (1) 

where q = p\ — p 2 , P = Pl + p 2 . The following notations are used: 7 s = i^ ^ 1 ^ 2 ^ 3 , 07^ = §[7^,7^], e 0123 = — 1, 
750" M y = — \tp,vapa a ^ ■, and S P ('f 5 j ,1 j L 'j a j^) = Aie' 1 " ' 13 . We study the form factors within the dispersion formulation 
of the quark model |l8). We start by considering the region q 2 < where the form factors may be represented as 
double spectral representations in the invariant masses of the initial and final qq pairs. The form factors corresponding 
to the decay region q 2 > are then derived by performing the analytical continuation. 

The transition of the initial meson q{rn 2 )q[m^) with the mass M\ to the final meson q(nii)q(m 3 ) with the mass M 2 
induced by the quark transition q(m 2 ) — > q(m±) through the current q(mi)0 fJl q(m 2 ) is described by the diagram of 
Fig.|l|. For constructing the double spectral representation we must consider a double-cut graph where all intermediate 
particles go on mass shell but the initial and final mesons have the off-shell momenta pi and and p 2 such that p\ = si 
and p\ = s 2 with (pi — p 2 ) 2 = q 2 kept fixed. 




FIG. 1. One-loop graph for a meson decay. 

The quark structure of the initial and final mesons is given in terms of the vertices Ti and r 2 , respectively. 
The initial pseudoscalar meson vertex has the spinorial structure T\ — 1756*1/ \/N c ; the final meson vertex has the 
structure T 2 — 175 G2/ ' \fN~c for a pseudoscalar state and the structure T 2ai = [A^^ + B{k\ — ks)^] G 2 /\fW c , A = — 1, 
B = 1/(t/s2~ + mi + 7713) for an S'-wave vector meson. 

The double spectral densities / of the form factors are obtained by calculating the relevant traces and isolating the 
Lorentz structures depending on pi and p 2 . These invariant factors / account for the two-particle singularities in the 
Feynman graph. 

For q 2 < the spectral representations of the form factors have the form 

00 sj(s 2 ,q ) _ ^ 

Z^ 2 ) = 7^ / ds 2^P2(s 2 ) I dsnpi(s x ) J 1 ) 2 S / 1 ' S2 ' q L , (2) 
16tH J J X i / 2 (si,s 2 ,q 2 ) 

(mi+m 3 ) 2 s~(s 2 ,q 2 ) 

where the wave function cpi(si) = Gj(sj) / (sj — M 2 ) and 

± 2 s 2 {m\ + m 2 — q 2 ) + q 2 {m\ + m%) — (m\ — m\){m\ — m 2 ) A 1 / 2 (s 2 , m§, m 2 )A 1 / 2 (g 2 , m 2 , m 2 ) 

and A(si, s 2 , S3) = (si + s 2 — S3) 2 — 4siS2 is the triangle function ^|. The analytical continuation of the expression (Q) 



2 The spectral densities / include proper subtraction terms. These subtraction terms have been determined in |18) by matching 
the structure of the heavy quark expansion in the quark model to the structure of the heavy-quark expansion in QCD 
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to the region q 2 > gives the form factor at q 2 < (ma — mi) 2 . Explicit expressions of the double spectral densities 
of all the form factors in ([j]) and more details can be found in [H [j 

The soft wave function of a meson M [q(m q )q(m q )] can be written as follows 



^s 2 - (m 2 q - m 2 ) 2 w(fc2) 
9?(S) ~ V2 J s -{m q -m q ) 2 ,3/4 W 

with k 2 = X(s,m 2 ,m^)/4:S. Here the ground-state radial S'-wave function w(k 2 ) is normalized as J w 2 (k 2 )k 2 dk = 1. 

As demonstrated in |L8j the form factors develop the correct structure of the long-distance corrections in accordance 
with QCD in the leading and next-to-leading 1/toq orders, if the radial wave functions w(k 2 ) are localized in a region 
of the order of the confinement scale, k 2 < A 2 . We assume a simple gaussian parametrization of the radial wave 
function 

w{k 2 ) oc exp(-fc 2 /2/3 2 ), (4) 

which satisfies the localization requirement for [3 ~ J^qcd- 

The leptonic decay constant of the pseudoscalar meson fp is given in terms of its wave function by the spectral 
representation Jl8| 

t /lv~ ( _i_ \ [ a I ^ 1/2 (s,m 2 q ,m 2 ) S - (m q - m q ) 2 

f P = ^/N c (m q + m q ) I ds <p(s) - . (5) 

/ an s s 



III. PARAMETERS OF THE MODEL AND NUMERICAL RESULTS 



A. Parameters of the model 

We consider the slope parameter ft of the meson wave function (Q) and the constituent quark masses as fit param- 
eters. The relevant values for the B, p, and it mesons have already been determined in |2Jj from a fit to the lattice 
results on the form factors T 2 (q 2 ) and ^4i(<? 2 ) (see Eq. (Q) below) at q 2 — 19.6 and 17.6 GeV 2 0. Thereby, use has 
been made of the spectral representation of the leptonic decay constant (|)|), and the double spectral representations 
(^J) of the form factors. The values obtained for m&, m u , and (3p, (3 P , Pw are displayed in Table ||. 

A few comments on these numbers are in order: 

• The quark model double spectral representations take into account long-range QCD effects but not the short- 
range perturbative corrections. However, by fitting the wave functions and masses to reproduce the lattice 
points, these corrections are effectively taken care of: Corrections to the quark propagators correspond to the 
appearance of the effective quark masses. Corrections to the vertices at the relevant values of the recoil variable 
uj = (M| + m 2 — q 2 )/2MpM^ should be small as found in form factors of other meson transitions [ p6j | . 

• The value obtained for the 6-quark mass nib — 4.85 GeV is close to the one-loop pole mass which in fact is the 
relevant mass for quark model calculations. 

We now need to fix the parameters describing the strange and charmed mesons. The charm quark mass can be 
determined from the well-known 1 /toq expansion of the heavy meson mass in terms of the heavy quark mass and the 
hadronic parameters A, Ai and A2. Using the recent estimates of these parameters one finds 

nib — Tn c ~ 3.4 GeV. (6) 

This provides m c ~ 1.45 GeV. For m s one expects m s ~ 350 — 370 MeV taking into account that m u = 230 MeV. 



3 The spectral representations (Q) take into account the long-distance contributions connected with the structure of initial and 
final mesons. To describe additional long distance effects, Eq. (^) should be multiplied by the constituent quark form factor 
/ 92 _> 91 (q 2 ) which contributes to the resonance structure in the q 2 channel. However, in the region of calculation q 2 < (mi — mi) 2 , 
the wave functions provide already for a rise of the form factors with q 2 , which is well compatible with a properly located meson 
pole. Thus, an additional quark form factor is not needed there, but we will use a proper extrapolation formula when considering 
the vicinity of the poles (see Eq (^) below). 
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A stringent way to constrain the parameters m s , (3k, (3k* , and (3d is provided by the measured integrated rates of 
the semileptonic decays D — ► (K,K*)lv. In addition we apply the relation (||) which connects (3k with m s by using 
the known value of the X-meson leptonic decay constant Jk — 160 MeV . The parameter values found this way are 
displayed in Table [jQ. The corresponding form factors and decay rates are given in Tables IV and [v| . 

The polarization of the K* in the D — > K*lv decays turns out to be in good agreement with the experimental 
result (Table VI), and the calculated D meson decay constant Jd = 200 MeV corresponds to the expectation for the 
magnitude of this quantity. 

The parameter (3d* cannot be found this way, but it should be close to (3d because of the heavy quark symmetry 
requirements. We therefore set (3d* = (3d- 

Also listed in Table | are the parameters which describe strange heavy mesons. They are discussed in subsection 

D. 



TABLE I. Constituent quark masses and slope parameters of the exponential wave function (in GeV). 



m u 


m a 


m c 


m b 


Ar 


(3 K 


(3d 


Pb 




Pd. 


Pb, 


(3 P 


Pk* 


(3d* 




0.23 


0.35 


1.45 


4.85 


0.36 


0.42 


0.46 


0.54 


0.45 


0.48 


0.56 


0.31 


0.42 


0.46 


0.45 



TABLE II. Leptonic decay constants of the pseudoscalar mesons in MeV calculated via H with the parameters of Table 1. 







fD 


fB 


fv. 


fn a 


fB. 


132 


160 


200 


180 


183 


220 


200 



TABLE III. Meson masses in GeV from PDG fl] 



M n 


M K 


M v 


M v , 


M D 


M Ds 


M B 


m Bb 


M p 


M K , 


M$ 


M D , 


M Dt 


0.14 


0.49 


0.547 


0.958 


1.87 


1.97 


5.27 


5.37 


0.77 


0.89 


1.02 


2.01 


2.11 



The knowledge of the wave functions and the quark masses allows the calculation of the form factors in Eq. (|lj) . It 
is however more convenient to present our results in terms of the dimensionless form factors F + , To, fx, V, Ao, A\, 
A 2 , T\ , T2 , T3 |l|] which are the following linear combinations of the form factors given in Eq. (|l|) : 

F+ = f+, F Q = f + + ^-f_, F T = -(M 1 + M 2 )s, 

V=(M 1 + M 2 )g, A x = M l + M f, A 2 =-(M 1 +M 2 )a+, A„ = (/ + q 2 ■ a_ + Pq ■ a+) , 

q 2 Pq 

Ti = -g+, T 2 = -g+-—g-, T 3 = <?_ - —g . (7) 

The form factors ^ are defined such that they involve only contributions of resonances in the q 2 channel with the 
same spin, whereas some of the form factors defined by the Eq. (Il|) contain contributions of resonances with different 
spins. The form factors F + , Ft, V, T\ contain a pole at q 2 = My = M 2 _ and Aq contains a pole at q 2 = Mb = M 2 _ 
(more details are given in the Appendix). 

The remaining form factors, Fq, A%, A 2 , T 2 and T3, do not contain contributions of the lowest lying negative parity 
states (for instance, F contains a contribution of the + state, and Ai contains that of the 1 + which have considerably 
higher masses). As a result they have a rather flat q 2 behaviour in the decay region, whereas the form factors F + , 
Ft, V, Ti, Aq are rising more steeply. 

From the spectral representations (|^) together with the parameter values of Table I the form factors are obtained 
numerically. For the applications it is convenient, however, to represent our results by simple fit formulas which 
interpolate these numerical values within a 1% accuracy for all q 2 values in the region < q 2 < (m 2 — mi) 2 . Also, 
they should be appropriate for a simple extrapolation to the resonance region. 



4 In J2f| a different set of the parameters was used which also provided a good description of the available experimental data 
on semileptonic B and D decays. However, the corresponding form factors have a rather flat q 2 -dependence and do not match 
the lattice results at large q 2 . 
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Let us start with the form factors F+, Ft, V, Ti, Ao. If we interpolate the results of the calculation with the simple 
three-parameter fit formula 



f(q 2 



/(0) 



{l-q 2 /M 2 )(l-q 2 /(aM) 2 )' 



(8) 



the least-% 2 interpolation procedure leads in all cases to a value of the parameter M which is within 3% equal to the 
lowest resonance mass. We consider this fact to be an important indication for the proper choice of the quark-model 
parameters and for the reliability of our calculations. We therefore prefer to fix the pole mass M to its physical value. 
The fit functions (||) represent the results now with an accuracy of less than 2%. To achieve the accuracy of less than 
1% in all cases we take the form [E4f: 



/(<z 2 ) 



/(0) 



(1 - q 2 /M 2 )[l - a iq 2 /M 2 + a 2 q 4 /M i ] ' 



(9) 



where M = Mp for the form factor Ao and M = My for the form factors F + , Ft, V, T\. In the Tables below we 
quote numerical values of a-i only if an accuracy of better than 1% cannot be achieved with 02 = 0> and take a% = 
if this accuracy can already be achieved with the two parameters /(0) and a\. A two-parameter fit was discussed in 

For the heavy-to-light meson transitions the masses of the lowest resonances are not very much different from the 
highest q 2 values in the decay. Eq. (||) then allows an estimate of the residues of these poles. These residues can be 
expressed in terms of products of weak and strong coupling constants (see Appendix) . The errors for these constants 
induced by changing o\ and 172 in our fitting procedure (keeping to the 1% requirement) do not exceed 10%. Moreover, 
the residues of the form factors at the meson pole are not independent and satisfy certain constraint (see Eq. ( [i"8| ) in 
the Appendix), which provides a consistency check of the extrapolations. The mismatch in ( |l8|) is always below 10%, 
and in most of the cases much lower. 

For the form factors Fq, A%, A2, T2 and T3 the contributing resonances (0 + , 1 + , etc) lie farther away from the 
physical decay region and the effect of any particular resonance is smeared out. For these form factors the interpolation 
formula taken is[] 

f(<l 2 ) = /(0)/[l - °iflMl + a 2 q 4 /Mfo. (10) 

If setting (72 = allows us to describe the calculation results with better than 1% accuracy for all q 2 , a simple 
monopole two-parameter formula is used. 

The values of /(0), ax, and <ii are given for each decay mode in the relevant subsections. 



B. Charmed meson decays 



1. D —> K, K* 



The D —f K, K* decays are induced by the charged current c — > s quark transition. As described in the previous 
section, the measured total rates of these decays are used for a precise fit of the parameters of our model. With the 
parameters of Table || we obtain the form factors listed in Table IV . Table |y| compares the form factors at q 2 = 
with the results of other approaches and Table VI presents the decay rates. 



TABLE IV. The D -» K,K* transition form factors. My = M D . = 2.11 GeV, M P = M Ds = 1.97 GeV. 
factors F+, Ft,V, Ao,Ti the fit formula Eq. (H) is used, for the other form factors - Eq. (|l(i|) 



For the form 





D -> K 


D -> K* | 




F+ 


Fq 


F t 


V 


Ao 


Ax 


A 2 


Ti 


T 2 


T 3 


/(O) 


0.78 


0.78 




0.75 


1.03 


0.76 


0.66 


0.49 




0.78 


0.78 


0.45 


171 


0.24 


0.38 




0.27 


0.27 


0.17 


0.30 


0.67 




0.25 


0.02 


1.23 


<72 




0.46 










0.20 


0.16 






1.80 


0.34 



5 One should note that the parameters ax and 02 in the fit formula jic| ) for the form factors Fo, Ax, A-x, T2, and T3 are 
introduced in a different way than in the fit formula (B|) for the form factors F+, Ft, V, Tx, and A$. 



TABLE V. Comparison of the results of different approaches for the semileptonic D — > K, K* form factors at q 2 = 0. 



Ref. 


F+(0) 


F T (Q) 


V(0) 


Mo) 


A a (0) 


This work 


0.78 


0.75 


1.03 


0.66 


0.49 


WQ"R |f{l 
VV D-D 


U. ( 




J..O 


U.oo 


1 9 


Jaus'96 |§ 


0.78 




1.04 


0.66 


0.43 


SR§ 


0.60(15) 




1.10(25) 


0.50(15) 


0.60(15) 


L at (average) [J13| 


0.73(7) 




1.2(2) 


0.70(7) 


0.6(1) 


Lat O 


0.71(3) 


0.66(5) 








Exp Q 


0.76(3) 




1.07(9) 


0.58(3) 


0.41(5) 



TABLE VI. The D -> (A", K*)lv decay rates in 10 1 



obtained within different approaches, |V CS | = 0.975. 



Ref. 


T(D — ► K) 


T(D — * A"*) 


r(A-*)/r(A-) 


r L /r T 


This v 
Jaus'9 
SR 


/ork 
3 § 


9.7 
9.6 
6.5(1.5) 


6.0 
5.5 
3.8(1.5) 


0.63 
0.57 
0.50(15) 


1.28 
1.33 
0.86(6) 


Exp 


3J 


9.3(4) 


5.7(7) 


0.61(7) 


1.23(13) Q 



Extrapolating the form factors to q 2 — M^, (or q 2 
constants (see the Appendix for the relevant formulas) 



9D'DKJv 



2M 



= 1.05 ±0.05, 



9d s dk»Jp 



M|, for Aq) gives the following estimates of the coupling 



D- 



2M 



K" 



= 1.7±0.1, 



AD*,) 

JT 

ADD 
Jv 



= 0.95 ±0.05, 



9d*dk' 
9d*dk 



= 1.1 ± 0.1. 



2. D 



7T, p 



These decays are induced by the c — > d charged current 
this mode allows for parameter-free predictions 



Since all the necessary parameters have already been f ixed, 
r II presents the results of our calculations. In Tables VIII 
and IX we compare our results with different approaches and with experimantal data. The form factors at q 2 = are 



close to the predictions of the relativistic quark model of Ref. m, but the q dependence is different such that our 
model and pj predict different decay rates. Although the experimental errors are very large and nearly all theoretical 
results agree with experiment, we notice perfect agreement of our decay rates with the central values. 



TABLE VII. The calculated D -> n,p transition form factors. M v = M D - = 2.01 GeV, M P = M D = 1.87 GeV . 
form factors F + , Ft, V, Aq,Ti the fit formula Eq. ^ is used, for the other form factors - Eq. (|Io|). 



For the 





D -> 7T 


I 




F+ 


F 


F T 


V 


A 


Ax 


A 2 


Ti 


T 2 




/(O) 


0.69 


0.69 


0.60 


0.90 


0.66 


0.59 


0.49 


0.66 


0.66 


0.31 




0.30 


0.54 


0.34 


0.46 


0.36 


0.50 


0.89 


0.44 


0.38 


1.10 


02 




0.32 














0.50 


0.17 



TABLE VIII. Comparison of the results of different approaches for the semileptonic D —* ir,p form factors at q 



Ref. 


F+(0) 


F T (0) 


V(0) 


Ax(0) 


A 2 (0) 


This work 


0.69 


0.60 


0.90 


0.59 


0.49 


WSB 


P 




0.69 




1.23 


0.78 


0.92 


Jaus'96 


\ 


0.67 




0.93 


0.58 


0.42 


SR 






0.50(15) 




1.0(2) 


0.5(2) 


0.4(2) 


Lat(ave) 


PI 


0.65(10) 




1.1(2) 


0.65(7) 


0.55(10) 


Lat | 


15 




0.64(5) 


0.60(7) 
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TABLE IX. The D -> (tt, p)lv decay rates in 10 10 s~\ \V cd \ = 0.22. 



Ref. 


F(D -» tt) 


r(F> -> P ) 


r( P )/i» 


r L /r T 


This work 


n op; 
U.yo 




U.40 


1 1 ^ 


WSB | 




0.68 


0.67 


1.0 


0.91 


Jaus'96 


k 


0.8 


0.33 


0.41 


1.22 


SR § 




0.39(8) 


0.12(4) 




1.31(11) 


Melikhov'97 |25[ 


0.62 


0.26 


0.41 


1.27 


Exp H 


0.92(45) 


0.45(22) 


0.50(35) 





For the coupling constants we get the following relations 

AD*) AD) ad*) 

9D " D ^v =1.05 ±0,05, 9DDpJp =2.1 ±0.2, q_=0.9±0.1, ^^=1.3 ±0.2. 

2M D , 2M P fiP*) g D , D7r 

Taking ff" ] ~ 220 MeV, we find 

gD'Dn = 18 ± 3 

in perfect agreement with a calculation of q d* d-k based on combining PC AC with the dispersion approach |3£ 



C. Beauty meson decays 



1. B ZX D* 



These decays arise from the heavy-quark b — > c transition. Here one has rigorous predictions for the expansion of the 
form factors in terms of the heavy-quark mass [ p0[ . Namely, the main part of the form factors can be expressed through 
the universal form factor - the Isgur-Wise function. However, different models provide different g 2 -dependences of the 
Isgur-Wise function as well as different subleading 1/itiq corrections. 

We recall that our spectral representations of the form factors explicitly respect the structure of the long-distance 
QCD corrections in the leading and the subleading orders of the heavy-quark expansion. Thus , we expect reliable 
predictions for the form factors. Our numerical results are summarized in Tables Kl XI, and XII. 



TABLE X. The calculated B — > D D* transition form factors. M v = M B * a 
F+,Ft, V, Ao, Ti the fit formula Eq. (H) is used, for the other form factors - Eq. 



Mp = M Bc = 6.4 GeV. 
&■ 



For the form factors 





B -> D 


II ■ 1)- | 




F+ 


F 


F T 


V 


Ao 


Ax 


A 2 


Ti 


T 2 


T 3 


/(0) 

02 


0.67 
0.57 


0.67 
0.78 


0.69 
0.56 


0.76 
0.57 


0.69 
0.58 


0.66 
0.78 


0.62 
1.40 
0.41 


0.68 
0.57 


0.68 
0.64 


0.33 
1.46 
0.50 



TABLE XL Comparison of the results of different approaches for the semileptonic B-tD, D* form factors at q 2 — 0. 



Ref. 


F+(0) 


Ai(0) 


Rx(Q) = V(Q)/Ax(0) 


R 2 (Q) = A 2 (0)/Ai(0) 


This work 
Jaus'96 Ml 
Neubert |p6| 
Close,Wambach p?[ 


0.67 
0.69 


0.66 
0.69 


1.15 
1.17 
1.3 
1.15 


0.94 
0.93 
0.8 
0.91 


Exp H 






1.18±0.15±0.16 


0.71±0.22±0.07 
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TABLE XII. The B (D, D*)lu decay rates in | V cb \ 2 ps' 1 . 



Ref. 


F{B -> D) 


T(B -> D*) 


r(U*)/T(£>) 


T L /r T 


This work 
Jaus'96 § 
Melikhov'97 25 1 


8.57 
9.6 
8.7 


22.82 
25.33 
21.0 


2.66 
2.64 
2.65 


l.n 

1.28 


Exp 


(1.34 ± 0.15) lO -2 ps~ L ||] 


(2.98 ± 0.17)10^ ps' 1 |l) 


2.35(1.3) 


1.24(0.16) |5]] 



For the coupling constants we find 
9B *c BDf v c) = i, 56 ±0,15, aB^BD-f^ =33±03 ^ = 0.9 ±0.1, ^££1 = 1.05 ±0.05. 



B-^K, K* 



These decays are induced by the 6 — ► s Flavor Changing Neutral Current (FCNC). We recall that the B — > 7r,p 
form factors at large g 2 have been used to fix the parameters of the model. Thus we expect that the predictions 
for the B — ^ K, K* form factors which in fact differ from the former mode only by SU(3) violating effects should be 
particularly reliable. Table XIII presents the calculated form factors and Fig ^ exhibits our predictions together with 
the available lattice results at large q 2 . The good agreement shows that the size and the sign of the SU(3) violating 
effects are correctly accounted for. 



TABLE XIII. The calculated B -> K,K* transition form factors. M v = M B * = 5.42 GeV, M P = Mb 
the form factors F+,Ft, V, Ao, Ti the fit formula Eq. (^) is used, for the other form factors - Eq. (|l"o|). 



5.37 GeV. For 





B -> K 


B -> K* I 




F+ 


F 


F T 


V 


Ao 


A 1 


A 2 


Ti 


T 2 


T 3 


/(0) 


0.36 


0.36 




0.35 


0.44 


0.45 


0.36 


0.32 




0.39 


0.39 


0.27 




0.43 


0.70 




0.43 


0.45 


0.46 


0.64 


1.23 




0.45 


0.72 


1.31 


02 




0.27 










0.36 


0.38 






0.62 


0.41 



TABLE XIV. Comparison of the results of different approaches for the B — * K, K* form factors at q 2 — 0. 



Ref. 


F+(0) 


F T (0) 


no) 


Mo) 


A 2 (0) 


Mo) 


21(0) 


r 3 (o) 


This work 


0.36 


0.35 


0.44 


0.36 


0.32 


0.45 


0.39 


0.27 


SR§_ 


0.25 




0.47 


0.37 


0.40 


0.30 


0.38 




Lat+StcchjH 






0.38 


0.28 




0.32 


0.32 




LCSR'98_ff 


0.34 


0.374 


0.46 


0.34 


0.28 


0.47 


0.38 


0.26 


Lat @| 


0.30(4) 


0.29(6) 















1.8 — m V 

: t a„ -■- 

1, 

u->*. i"; 

o.s — 

o.4 : ;;;'. ;; _ v ., -••*♦•• 

0.2 — ^-^^ ? tk~ 







10 



15 



q 



FIG. 2. Form factors of the B — > K* transition vs the lattice results. 
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For the coupling constants we obtain 

9B ''** f v s) = 0-65 ±0.05, 9B '**'/* m) = 1-65 ±0.1, 4tJt = °-95 ±0.05 = 1.15 ± 0.05. 

2M B; 2Mk" ffs) g B . BK 



3. B — > 7T, p 



The B —> p transition has been used for determining the parameters of our quark model in the u, d and b sectors by 



fitting the quark-model form factors to available lattice results on T 2 and A\ at large q in 



The corresponding 



form factors and the decay rates have been calculated in this art icle. We present the results of |24] in terms of 
parametrizations for the form factors of the set (|t]) (see Table |XV| ) which have not been given in that paper. The 
only difference of the results presented here with the results of p4| occurs for the B — > tt form factors F + and Fq at 
q 2 > 20 GeV 2 . In H] these quantities are not calculated directly but extrapolated from the region q 2 < 20 GeV 2 . 
The parametrizations of /+ and /_ in p4j correspond to gB*Bw — 50. The parametrization given here, on the other 

Ai • 20 GeV 2 



hand, corresponds to gB*Bn — 32 which is in agreement with recent analyses of this quantity 
both parametrizations describe the results of the numerical calculation well and agree with the available results from 
lattice QCD at q 2 < 22 GeV 2 Q. For these reasons, the earlier result T(B — > irlv) — 8.3|Vu&| 2 ps 1 calculated with 



the parametrization of the form factor F+ from Table XV remains practically unchanged compared to 

The form factor Fq at large q 2 lies below the central lattice values but nevertheless agrees with lattice results within 
the given error bars. Notice however that in our model the form factor F is calculated as a difference of / + and /_ 
and at large q 2 turns out to be much more sensitive to the subtle details of the pion wave function, than / + and /_ 



separately. A simple Gaussian wave function which works quite well for / + and /_ 
for F . 



might not be sufficiently accurate 



TABLE XV. The calculated B — > n,p transition form factors. M v = Mb* = 5.32 GeV, M P = M E 
form factors F + , Ft, V, Aq, T\ the fit formula Eq. (§) is used, for the other form factors - Eq. (|lo|). 



= 5.27 GeV. For the 





B — > 7T 


B^p | 




F+ 


F 


F T 


V 


A 


Ai 


A 2 


Ti 


T 2 




/(0) 


0.29 


0.29 


0.28 


0.31 


0.30 


0.26 


0.24 


0.27 


0.27 


0.19 




0.48 


0.76 


0.48 


0.59 


0.54 


0.73 


1.40 


0.60 


0.74 


1.42 


02 




0.28 








0.10 


0.50 




0.19 


0.51 



Table [XVI compares the results obtained from the quark model of Ref. 0| with results from the quark model 
of Jaus Q and latest light-cone sum rule (LCSR) results pi] . One observes very good agreement between the 
quark model of Jaus, LCSRs, and our approach. The only visible difference with the LCSR method occurs in 
the form factor Ao(0), which is caused by small differences of the two methods in Ai(0) and ^(0) (recall that 
Aq(0) = ((Mi +M 2 )Ai(0) - (Mi - M 2 )A 2 (0)) /2M 2 ) This discrepancy exceeds the 15% error bar quoted for the 
LCSR results only marginally. If the LCSR results at small q 2 and lattice results at large q 2 are correct, our approach 
surely provides a realistic description of the form factors at all kinematically accessible q 2 values. 



TABLE XVI. Comparison of the results of different approaches on weak B — > tt, p form factors at q 2 — 0. 



Ref. 


F+(0) 


F T (0) 


V(0) 


Mo) 


Mo) 


Mo) 


Ti(0) 


n(o) 


This work 


0.29 


0.28 


0.31 


0.26 


0.24 


0.29 




0.27 


0.19 


Jaus'96 fll 


0.27 




0.35 


0.26 


0.24 










LCSR'98jl) 


0.305 


0.296 


0.34 


0.26 


0.22 


0.37 




0.29 


0.20 


Lat ||| 


0.28(4) 


0.28(7) 

















Extrapolating the form factors to the poles, we obtain 



AB*) 

gB-BnJy 

2M B ~ 



= 0.6 ±0.05, 



2M n 



(B) 

p =1.4 ±0.2, 



Using fy ~ 200 MeV gives the estimate 



gBB'-K — 32 ± 5, 



f(B') 
JT 
AB') 
JV 



= 0.97 ±0.03, 



SB' Bp 
gB*Bn 



= 1.2 ±0.1. 



2y/M B M B * 



g B B* w = 0.4 ±0.06 



(11) 
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The latter value is in good agreement with the lattice result g = 0.42 ± 0.04 ± 0.08 ^d[] and is only slightly smaller 
than g = 0.5 ± 0.02 |3j| based on combining PCAC with our dispersion approach. 

Summing up our results on the decays of the nonstrange heavy mesons, we found no disagreement neither with the 
exisiting experimental data nor with the available results of the lattice QCD or sum rules in their specific regions of 
validity. The only exception is the form factor Fq at large q 2 in B — > ir and D — > ir decays, where our results are 
lying slightly below the lattice points. However, this disagreement can be related to a strong sensitivity of Fq at large 
q 2 to the details of the pion wave function. Small changes in the pion wave function, which only marginally affect /+ 
and /_, can change the form factor Fq. But such subtle effects are beyond the scope of our present analysis. 

In the next section we apply our model to the decays of strange heavy mesons for which a few new parameters have 
to be introduced which are specific to the description of weak decays of strange heavy mesons to light mesons. 



D. Decays of the strange mesons D s and B s 



Before dealing with these decays, we must first specify the slope parameters of the B s and the D s wave function. 
We obtain these parameters by applying (^|) and using /b,//b = 1.1 and /d s //d = 1.1 in agreement with the lattice 
estimates for these quantities jl3| . The resulting values of the slope parameters are listed in Table ffl. Since all other 
parameters have already been fixed the calculation of the form factors is straight forward. The only exceptions are 
the decays into the 77, rj', (j) final states. For these decays we need to know the <f> wave function, the mixing angle and 
the slope of the radial wave function of the ss component in 77 and 77'. Our procedure of fixing these parameters are 
discussed in the relevant subsection. 



1. D a 



K, K* 



These meson transition are driven by the charged-current 



d quark transition. The results of the calculation 



are given in Table XVII. The predictions for the semileptonic decay rates are displayed in Table XVIII 



TABLE XVII. The calculated D a -» K, K* transition form factors. M v = M D * = 2.01 GeV, M P = 
the form factors F+, Ft, V, Aq,Ti the fit formula Eq. (^) is used, for the other form factors - Eq. (|lo|). 



M D = 1.87 GeV. For 





D S ^K 


/ >. • l< ■ I 




F+ 


Fo 


Ft 


V 


A 


Ai 


A 2 


Ti 


T 2 


Ta 


/(0) 


0.72 


0.72 




0.77 


1.04 


0.67 


0.57 


0.42 




0.71 


0.71 


0.45 




0.20 


0.41 




0.24 


0.24 


0.20 


0.29 


0.58 




0.22 


-0.06 


1.08 


02 




0.70 










0.42 








0.44 


0.68 



TABLE XVIII. The D s -> (K, K*)lv decay rates in 10 10 a -1 , \V cd \ = 0.22. 



Ref 


V{D S -> K) 


T{D S -» K*) 


T{K*)/T{K) 


IV/Fr 


This work 


0.63 


0.38 


0.6 


1.21 



For the coupling constants we obtain 

9 D ' D - K fv I =0.95 ±0.05, 9DD.K*f { f> =L85±(U5 4^= 0.9 ±0.1, = 1.15 ± 0.15. 

2M D , 2M K , f(P') g D *D s K 
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2. D.^ri,r)',4> 



These decay modes are induced by the charged current c — > s quark transition. The pseudoscalar mesons 77 and 77' 
are mixtures of the nonstrange and the strange components with the flavour wave functions r\ n = uu +^ d and rj s — ss, 
respectively, 



with the angle <p ~ 40° 



r) = coa(ip) r] n - sm(ip) r) s 
i]' = sm((p) rjn + cos((p) r) s , 

The decay rates of interest are 

T(D S -> r]lv) = sin 2 (<p)T(D a -> Vs(M n )lv) 
T{D S ri'lu) = cos 2 {v)T(D s r] s {M T} i)lv). 



(12) 



(13) 



Let us give a brief explanation of these formulas: The semileptonic decay rates are determined by the form factor / + . 
The spectral representation of this form factor does not involve the final meson mass explicitly. This means that for 
the ss component of both 77 and 77' we have to deal with the same form factor, which can be expressed through the 
radial wave function of this component. On the other hand, the phase-space volume of the decay process is determined 
by the physical meson masses, as indicated in (p"3|). It should be clear, however, that the rj s is not an eigenstate of 
the Hamiltonian and does not have a definite mass. 

Assuming universality of the wave functions of the ground state pseudoscalar 0~ nonet, the radial wave function 
of the nonstrange component ^-q n coincides with the pion radial wave function f42|| . The radial wave function 
should be determined independently. From the analysis of a broad set of processes the leptonic decay constant f a of 
the strange component rj s , has been found to lie in the interval f s = (1.36± 0.04)/^ [jt3[ . This allows us to determine 
the slope parameter (3 rja in such a way that the calculated value of f s lies in this interval, and the calculated ratio 
T(D S — > rj)/T{D s — > 7/) agrees with the experimental data for ip = 40°. This procedure yields for the slope parameter 
(3 rjs = 0.45 M. For the slope parameter of the wave function of the </>-meson, which is the vector ss state, we expect 
a value close to fl Vs . 

In fact, /?0 = 0.45 GeV leads to the B s — > (j> transition form factors which agree well with the LCSR results at 
q 2 — (see subsection 4). With all other quark model parameters fixed from the description of the nonstrange heavy 
meson decays and by taking a simple Gaussian form of the radial wave function, the decay rate T(D S — > 4>lv) is a 
function of (3$. This function has a minimum at the value /3$ — 0.45 GeV; nevertheless, the corresponding value of 
the decay rate is la above the central experimental value ) . 

The results of our calculations are given in Tables [XIX| and |XX[ . 



TABLE XIX. The calculated D s -» 77^ transition form factors. M v = M D j = 2.11 GeV, M P = Md, = 1.97 GeV. For 
the form factors F+,Ft, V, Ao,Ti the fit formula Eq. (Jojj is used, for the other form factors - Eq. (|Io|). 







D s -> Vs(M v ) 


D a -> T) s (M v i) 


IK. > I 




F+ 


Fa 


F T 


F 


F T 


V 


A 


A 1 


A 2 


Ti 


T 2 


T 3 


/(O) 


0.78 


0.78 


0.80 


0.78 


0.94 


1.10 


0.73 


0.64 


0.47 


0.77 


0.77 


0.46 




0.23 


0.33 


0.24 


0.21 


0.24 


0.26 


0.10 


0.29 


0.63 


0.25 


0.02 


1.34 


02 




0.38 




0.76 














2.01 


0.45 



TABLE XX. The D s -» (77, 77', <f)lv decay rates in 10 10 s" 1 , |V CS | = 0.975. The experimental rates are obtained from the 
corresponding branching ratios using the D s lifetime td 3 = 0.495 ± 0.013 ps from the 1999 update |3l| 



Ref 


T(D S -» 7 ? ) 


F(D 3 -» rf) 


T(D S c/>) 


This work 
Exp H 


5.0 
5.2±1.3 


1.85 
2.0±0.8 


5.1 
4.04±1.01 



3 Another procedure of taking into account the SU(3) breaking effects to obtain \[' r|s from \£ r r)n has been proposed in 
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For the coupling constants we obtain 
9DiDsvJv s =L0±(U 9D ° D J?! ipDs) =1-6 ±0.3, 4S! = 0.93 ±0.03, = 1.08 ± 0.04. 

2M]j* 2M <j) f(p:) g D * DsVs 



3. B s -> K, K* 



This mode is driven by the b —> u charged current transition. The only additional new parameter needed here is 
the slope of the B s wav e function. We obtain it by using (||) and taking JbJJb = 1.1. The results of our calculation 
are given in Table XXI. 



TABLE XXI. The calculated B s -» K, K* transition form factors. M v = M B * = 5.32 GeV, M P = M L 
the form factors F+,Ft, V, Aq,T± the fit formula Eq. (^) is used, for the other form factors - Eq. (|i"o|). 



5.27 GeV. For 











F+ 


F 


F T 


V 


A 


A x 


A 2 


7i 


T 2 


T 3 


/(0) 


0.31 


0.31 




0.31 


0.38 


0.37 


0.29 


0.26 




0.32 


0.32 


0.23 




0.63 


0.93 




0.61 


0.66 


0.60 


0.86 


1.32 




0.66 


0.98 


1.42 


(72 


0.33 


0.70 




0.30 


0.30 


0.16 


0.60 


0.54 




0.31 


0.90 


0.62 



These form factors lead to the following relations 



,(B*) JB) 

9b* b„k Jy _ =0<44±0>04i 9BB - K ; f r =1.3 ±0.1, 



2M T _ 



2M 



K* 



f(B*) 

JT 

AB') 

Jv 



0.95 ±0.05, 



9b*b 3 K' 
9b'B b k 



= 1.2 ±0.1. 



The form factors at q = are compared with the LCSR predictions in Table XXII. We observe some disagreement 



between our predictions and the LCSR calculation which gives smaller values for all the form factors. A closer look 
at the origin of this discrepancy shows that its source is the strength and sign of the SU(3)-breaking effects. They 
lead to opposite corrections in the two approaches. 



TABLE XXII. Comparison of the QM and LCSR results on the B s -> K, K* form factors at q 2 = 0. 



Ref. 


V(0) 


Mo) 


A 2 (0) 


A o (0) 


Ti(0) 


T 3 (0) 


This work 
LCSR'98 


0.38 
0.262 


0.29 
0.19 


0.26 
0.164 


0.37 
0.254 


0.32 
0.22 


0.23 
0.16 



To discuss these SU(3) breaking effects, let us start with B — > p, which in fact differs from the B s — > K* only by 
the flavour of the spectator quark, and move to B s — > K* by accounting for the SU(3) violating effects: 

Within the LCSR method there are two changes which affect the form factors: first, the change Jb s — > Jb leads 
to an increase of the B s — » K* form factors; second, the change of the symmetric twist-two distribution amplitude of 
the /o-meson to the asymmetric one of the K* meson leads to a decrease of the form factors. The second effect turns 
out to be much stronger than the first one with the result of an overall decrease of the form factors. 

In the quark model, the same SU(3) breaking effects take place: The change of the spectator mass (it determines 
the increase of /b s //s) and the change of the K* meson wave function (due to the change of both the quark mass 
and the slope parameter of the light meson wave function). Here the influence of the slope of the heavy meson wave 
function upon the form factor is only marginal. Therefore, the resulting effect of these changes leads to an increase 
of the form factors. 

We want to point out that the difference between the results of the two approaches does not arise from specific 
effects (higher twists, higher radiative corrections etc) which are present in the LCSRs but absent in the quark model. 
The observed difference is only due to the different strength of the SU(3) violating effects at the level of the twist-2 
distribution amplitude. As was discussed in |26[ , this distribution amplitude can be expressed through the radial soft 
wave function of the meson. The change of the quark-model wave function caused by SU(3) violating effects does not 
induce a strong asymmetry in the leading twist-2 distribution amplitude. 

In view of the discrepancy between our results and the LCSR it would be interesting to have independent calculations 
of the B s — » K* form factors at small q 2 from the 3-point sum rules, as well as a lattice calculation for large q 2 . 
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4- B a -> r},r]',(j> 



The results of th e form factor 
2 _ n T n w n RFytvI The 



These weak meson transitio ns are i nduced by the FCNC b — > s quark transition, 
calculation are given in Table XXIII and compared with the LCSR predictions at q 2 = in Table XXIV 
agreement between the two values is satisfactory at least within the declared 15% accuracy of the LCSR predictions. 
This allows us to expect that also the D s — > <fi, r/, r}' form factors and the corresponding decay rates (given earlier in 
subsection 2) are calculated reliably. 



TABLE XXIII. The calculated B s -> 77^, transition form factors. M v = Mg* = 5.42 GeV, M P = M Bs = 5.37 GeV. For 
the form factors F+ , Ft , V, Ao , Ti the fit formula Eq. (^) is used, for the other form factors - Eq. (|lo|). 







B a - 


-> rj s (M v ) 


B s - 




B s ^<j> I 




F+ 


F 


Ft 


Fa 


F T 


V 


Ao 


Ai 


A 2 


Ti 


T 2 


T 3 


/(0) 


0.36 


0.36 


0.36 


0.36 




0.39 


0.44 


0.42 


0.34 


0.31 


0.38 


0.38 


0.26 




0.60 


0.80 


0.58 


0.80 




0.58 


0.62 


0.55 


0.73 


1.30 


0.62 


0.83 


1.41 


02 


0.20 


0.40 


0.18 


0.45 




0.18 


0.20 


0.12 


0.42 


0.52 


0.20 


0.71 


0.57 



TABLE XXIV. Comparison of the QM and LCSR results on the B B — > (f> form factors at q 2 = 0. 



Ref. 


V(0) 


Mo) 


A 2 (0) 


A o (0) 


Ti(0) 


r 3 (o) 


This work 


0.44 


0.34 


0.31 


0.42 


0.38 


0.26 


LCSR'98 


0.433 


0.296 


0.255 


0.382 


0.35 


0.25 



For the coupling constants we obtain 

;b sV J v =06±0Q5 9bj _ _ 

2M B: ' 2M f (Bt) ' 9 b* b 



9b:b.vJv' =0 .6±0.05, 9B 3 B a 4,f P Bs) =L5±QJ = 0.95 ± 0.05, = 1.13 ± 0.06. 



IV. CONCLUSION 



We have calculated numerous form factors of heavy meson transitions to light mesons which are relevant for the 
semileptonic (charged current) and penguin (flavor-changing neutral current) decay processes. Our approach is based 
on evaluating the triangular decay graph within a relativistic quark model which has the correct analytical properties 
and satisfies all known general requirements of QCD. 

The model connects different decay channels in a unique way and gives the form factors for all relevant q 2 values. 
The disadvantage of the constituent quark model connected with its dependence on ill-defined parameters such as the 
effective constituent quark masses, have been reduced by using several constraints: the quark masses and the slope 
parameters of the wave functions are chosen such that the calculated form factors reproduce the lattice results for the 
B — > p form factors at large q 2 and the observed integrated rates of the semileptonic D — > K, K* decays. 

Our main results are as follows: 

• In spite of the rather different masses and properties of mesons involved in weak transitions, all existing data 
on the form factors, both from theory and experiment, can be understood in our quark picture. Namely, all 
the form factors are essentially described by the few degrees of freedom of constituent quarks, i.e. their wave 
functions and their effective masses. Details of the soft wave functions are not crucial; only the spatial extention 
of these wave functions of order of the confinement scale is important. In other words, only the meson radii are 
essential. 

• The calculated transition form factors are in all cases in good agreement with the results available from lattice 
QCD and from sum rules in their specific regions of validity. The only exception is a disagreement with the 
LCSR results for the B s — > K* transition where we predict larger form factors. This disagreement is caused 
by a different way of taking into account the SU(3) violating effects when going from B — * p to B s — > K* 
and is not related to specific details of the dynamics of the decay process. We suspect that the LCSR method 
overestimates the SU(3) breaking in the long-distance region but this problem deserves further clarification. 
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• We have estimated the products of the meson weak and strong coupling constants by using the fit formulas for 
the form factors for the extrapolation to the meson pole. The error of such estimates connected with the errors 
in the extrapolation procedure is found to be around 5-10%. 

We cannot provide for definite error estimates of our predictions for the form factors because of the approximate 
character of the constituent quark model. However from the fine agreement obtained in cases where checks are 
possible, we believe that the actual accuracy of our predictions for the form factors is around 10%. Since some 
parameters have been fixed by using lattice results and have also been tested using the sum rule predictions, further 
improvements of the accuracy of our predictions will follow if these approaches attain smaller errors. Of course, each 
precisely measured decay will also allow a more accurate determination of the parameters of the model and thus can 
be used to diminish the errors at least for closely related decays. 
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VI. APPENDIX: WEAK AND STRONG MESON COUPLING CONSTANTS 

We provide here definitions of the coupling constants which determine the behavior of the form factors at q 2 near 
the resonance poles (beyond the decay region). We consider as an example the case of the B — > ir, p transition. 

1. Weak decay constants 

Weak decay constants of mesons are determined by the following relations 

(B( g )|6(0) 7M 759(0)|0)=zjfV 

(B% q )\HQ) 1 , q (m=t { r ) M B *f { v B ' ) 
(B*(q)\b(0)a^q(0)\0) = i(e^q v - e^q^)f T B '\ (14) 



where ^ is the B* polarization vector. In the heavy quark limit one has f p B ^ = /y B ^ = f T B ^ 

2. Strong coupling constants 
Strong coupling constants are connected with the three-meson amplitudes as follows 

(n(p 2 )B( Pl )\B*(q)) = \g B , B ^f) 
(p(p 2 )B( Pl )\B*(q)) = l -e a ^e { pP, qv 9 -^ 

(p(p 2 )B( Pl )\B*(q)) = l -g B *B P P^\ (15) 

where q — pi — p 2 , P = Pi + Vi an d is the polarization vector of the vector meson. 
In the heavy quark limit g B * Bp = 9bb p - 



3. Form factors 

The form factors F + , Ft, V, T\ contain pole at q 2 — M B , due to the contribution of the intermediate B* (1~ 
state) in the q 2 channel. The residue of this pole is given in terms of the product of the weak and strong coupling 
constants such that in the region q 2 ~ M B « the form factors read 
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t(B') 

9B'B7rJy 



2M* l-q 2 /M 2 



_ 9B'B,! { T B ' ] Mb + Mk 1 
T 2M* M* B l-q 2 /M 2 , ' 

9B*B P f ( v B * ] M B + M p 1 

2M* M B . l-q 2 /M%, "" 

T = 9B'B P Jt 1 , 

1 2M* l-q 2 /M%, 

(16) 

Here . . . stand for the terms non-singular at q 2 = Mj>* . 

Similarly, Aq contains the contribution of the B (0~ state). In the region of q 2 ~ M B it can be represented as 
follows 

, 9BB P f ( P B) M B 1 
°~ 2M B 2M p l-q 2 /M 2+ - { > 

First let us notice that the residues of the form factors are not all independent and are connected with each other as 
follows: 



Res{F T )Res(V) _ M B + M p M B + M w 
i?es(F + ) J Res(Ti) ~~ ~M B ~, M^ 



(18) 



This relation can be used as a cross-check of the consistency of the extrapolation for the form factors. 
The coupling constants are related to the residues of the form factors according to the relations 

( (B*) ,(B) 
9B*BttJ v d ,„ 9BBpJp nr> x 

2Mb! = ReS{F+) ' = 2i?6S(Ao) ' 

f [ T B ' ] _ Res(F T ) M B , g B , Bp _ Res(V) M B * 



fCB') Res(F+) M B + g B , B7T Res{F+) M B + M p 



(19) 
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